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1 Introduction 



The production of high-p-^ jets probes the short-distance dynamics of photon-photon reactions. 
In addition to providing tests of perturbative QCD, 77 processes with real or almost real 
photons give us information on the photon structure function which is complementary to the 
information gained from deep inelastic scattering on a real photon. The latter process essen- 
tially probes the quark distribution while high-p T jet production is also sensitive to the gluon 
distribution of the photon. 

As it is well known high-p^ processes induced by real photons have a rather complex struc- 
ture. This arises from the fact that the photon couples to the hard subprocess either directly 
or through its quark and gluon content. 

In leading order QCD (LO) the cross section for the production of jets can therefore be 
decomposed into a direct, single- resolved, and double-resolved component. In next-to-leading 
order (NLO) the photon-quark collinear singularities in the direct and single resolved compo- 
nents are subtracted and absorbed into the quark distribution of the photon in accord with 
the factorization theorem. This subtraction at the factorization scale M introduces an interde- 
pendence of the three components so that the separation into the direct, single-resolved, and 
double-resolved contributions becomes scale dependent. This means that in NLO all three com- 
ponents must be considered together and consistently be calculated up to NLO in the photon 
structure functions and in the hard scattering cross sections. 

Complete NLO calculations have been done previously for the case of the inclusive single-jet 
cross section [|I], ^ |3j and compared to experimental data from TRISTAN ||, 0]. Satisfactory 
agreement between data and theory has been found. In our previous work 0, |3[] we also calcu- 
lated the NLO inclusive two-jet cross sections for the direct and single-resolved contribution and 
estimated the double-resolved contribution by the LO calculation with k factors taken from the 
NLO inclusive one-jet cross section. In the meantime the calculation of the NLO corrections for 
the double-resolved cross section has been completed also for the inclusive two-jet cross section 
using techniques also employed for the direct and single-resolved cross sections in our previous 
work[|,|. ' 

In this paper we shall present the complete NLO predictions for the inclusive two-jet cross 
section incorporating all three components, direct, single-resolved and double-resolved. We 
shall compare these predictions with the older TRISTAN data |4], || and with recent data from 
LEP at e+e" energies of 130 and 136 GeV (LEP1.5) [§. The details of the NLO calculation of 
the double- resolved cross section are given elsewhere 0. Those of the direct and single- resolved 
cross section are described in 

The outline of the rest of the paper is as follows. In section 2 we explain the formalism 
used to calculate the inclusive two-jet cross section. Section 3 contains our results for the 
cross sections and some studies on the overall scale dependence and the dependence of the 
two-jet cross section on the cone radius R. Here we also present the results for different sets of 
photon structure functions and the comparison with the TRISTAN and LEP1.5 data. Section 
4 contains a short summary. 

2 Inclusive Two-Jet Cross Sections 

As explained in the introduction inclusive jet cross sections at large pt receive contributions 
from three parts according to how the incoming photons take part in the hard scattering 
subprocesses. For example, when both photons are resolved the inclusive two-jet cross section 
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has the following form 

d 3 a DR 
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In (|T|), pt is the transverse momentum of the measured or trigger jet with rapidity 771. 772 
denotes the rapidity of another jet such that in the three-jet sample (in the NLO calculation 
only up to three jets occur in the final state, the photon remnant jets are not counted) these two 
jets have the largest and second largest pr (pti,Pt2 > Pt$)- The first term on the right-hand 
side of (HD stands for the LO cross section and the second term is the NLO correction which 
depends on the factorization scale M at which the initial state singularities are absorbed into 
the parton distributions of the photon. /1 is the renormalization scale. The variable R is the 
usual cone size parameter which defines the size of the jets with transverse momentum p?, 
rapidity rj and azimuthal angle <fi. When two partons fulfill the Snowmass constraint || with 
the cone size parameter R they are recombined in one jet. The same jet definition is used in the 
analysis of the experimental data by the TOPAZ || , AMY |J and OPAL M collaborations. It 
is clear that either of the two jets can consist of a recombination of two partons inside a cone 
with radius R. 

The NLO corrections in the hard scattering processes are calculated with the phase space 
slicing method using for the separation of the 77 — > 2 jets and 77 — > 3 jets cross sections an 
invariant mass cut-off y, defined as 2piPj < ys, where s is the partonic center-of-mass energy 
squared. For example, the cross section for the direct process 77 — > qqg has soft, initial, and 
final state collinear singularities. The 77 — > qqg cross section is integrated over these singular 
regions up to the cut-off y which isolates the respective singularities. These cancel against the 
singular contributions of the virtual corrections to 77 — * qq and the subtraction terms at the 
scale M to be absorbed into the parton distribution functions of the photon. Outside the cut-off 
region controlled by the parameter y we have genuine qqg final states. In this contribution two 
of the partons are combined if they obey the Snowmass constraint. The LO contribution, the 
NLO virtual contributions, and the NLO corrections inside the y cut contribute to the two-jet 
cross section together with the contributions inside the cone with radius R. The part of the 
qqg cross section not fulfilling the cone recombination condition is the 3-jet cross section from 
which we have calculated the inclusive 2-jet cross section as a function of pr, T]i and rj2. 

For the exclusive 2-jet cross section p Tl = p T2 = p T . To separate double singularities in the 
qqg contributions we have used the method of partial fractioning. The same partial fractioned 
expressions are applied for the calculation inside and outside the y cut. Inside the y cut the 
calculation is done analytically in the approximation that terms 0(y) are neglected. Outside the 
y cut the contributions to the inclusive cross section are evaluated numerically with no further 
approximations. Because of the approximation in the analytic part, the parameter y must be 
chosen very small, of the order of 10~ 3 to 10 -4 . For such values of y the inclusive cross section 
should be independent of the invariant mass parameter y. However the inclusive two-jet cross 
section depends on the cone size R, which must be chosen as in the analysis of the experimental 
data. In the case of the single- and double-resolved contributions the calculation of the NLO 
corrections proceeds in the same way. For example, in the single-resolved case, we have in LO 
the contributions from jg — > qq and jq — > qg and in NLO we must consider jg — > qqg, jq — > qqq 
and 75 — > qgg which are treated in the same fashion as the 3-parton contributions in the direct 
case. In the calculation of the double-resolved cross section we have no photons in the initial 
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state, only quarks and gluons, and we have many more channels to consider as processes with 
three partons in the final state, which will not be given here. Otherwise the calculation proceeds 
in a completely analogous way to the direct and single-resolved case. 

Before the final results were obtained some tests of the NLO corrections have been per- 
formed. First we checked that the cross sections are independent of the cut-off y if y is chosen 
small enough. This was the case for y < in all considered cases. For y > 10~ 3 we observed 
some small y dependence which is caused by the approximation in the analytical part. Sec- 
ond the inclusive one-jet cross section for the double-resolved contribution was calculated and 
compared with our earlier results for which a completely different method for the cancellation 
of the infrared and collinear singularities was applied 0, |3|] . Very good agreement was found in 
all channels. Third we tested that the sum of the NLO direct and the LO single-resolved cross 
section is independent of the factorization scale M. The same test was performed for the sum 
of the NLO single-resolved cross section and the LO double-resolved cross section. Details of 
these tests can be found in 0. In this work also the details of the analytic calculation and the 
cancellation of the singularities for the NLO corrections of the double-resolved cross section 
are given. For the direct and the single-resolved cross sections the calculation of the NLO 
corrections was described earlier in [Q and 0, respectively, and also in 0. 

For the calculation of the single- and double-resolved cross sections we need the parton 
distributions of the photon ify 7 which appear in (p]) and in the equivalent formula for the 
single-resolved cross section (see ||). We have chosen the NLO set of Gliick, Reya and Vogt 
(GRV) in the MS scheme |TIJ as our standard set. For a comparison we shall employ also 



the MS sets of Gordon and Storrow |LT| and of Aurenche et al. [Q. We choose all scales 
H = M = p T except when we test the scale dependence. a s is calculated from the two-loop 
formula with N f = 4 (5) massless flavours with = 0.20 GeV(A^ = 0.13 GeV) for the e+e" 
center of mass energy \fS = 58 GeV(\/~S = 133 GeV). The charm, respectively, bottom quarks 
are also treated as light flavours with the boundary condition that the charm (bottom) content 
of the photon structure function vanishes for M 2 < m 2 {m 2 ) with m c = 1.5 GeV(m& = 5.0 GeV) 
for the GRV parton distributions [|n| . For the other two sets the treatment of heavy flavours 
is somewhat different. The details are found in the respective references [[TT], [T2] ]. 



3 Results and Comparison with Data from TOPAZ, AMY, 
and OPAL 

In e + e~ collisions the jets are produced via the exchange of two almost real photons. The 
spectrum of these small Q 2 photons is described by the Weizsacker- Williams approximation 



7 a l + {l-z) 2 f J E 2 9 2 max (l - z) 2 + my 
= 27 z ln ml? 



e 



m 2 z 



+ 2(1 Z) [E 2 e 2 max (l-z) 2 + m 2 e z 2 2)) (2) 

where ^ max is the maximally allowed angle of the electron (positron) which is given experi- 
mentally by the anti-tagging counters. E is the beam energy of the electron (positron) and 
z = E^/E is the fraction of the electron (positron) energy transferred to the respective photon. 
Then the inclusive two-jet cross section for the reaction e + + e~ — > e + + e~ + jet x + jet 2 + X 
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is calculated from 

d 3 cr(e + + e~ -> e + + e~ + jet x + jet 2 + X) 
dp T d?7id?72 

/• dV( 7 + 7 -> jet! + jet 2 + X) 

dzi dz 2 F l/e {zi,E)F 1 / e {z 2 ,E) — - — (3) 

J dp T dr]idri2 

where the cross section on the right-hand side is for the interaction of two real photons (Q 2 = 0) 
and is obtained from ([!]) for the double-resolved case and similar formulas for the direct and 
single-resolved 77 cross section (see ||). 

First we show in Fig. 1 the dependence of the inclusive 2-jet cross section on the cone radius 
R for the special case of LEP2 kinematics, y/S = 175 GeV, # max = 1.72°, since we cannot expect 
this dependence to be tested with TRISTAN or LEP1.5 data. Of course, the LO cross section is 
independent of R which is shown as the broken line in Fig. 1. In the NLO calculation the cross 
section increases with increasing R, first almost linearly in InR for R values up to 0.75 and 
somewhat stronger for R > 0.75. At R ~ 0.9 the LO and NLO cross sections are equal for the 
special kinematical point px = 10 GeV, 771 = 772 = 0. Therefore for R values, chosen between 
0.7 and 1.0, the higher order corrections are moderate. The LO curve is obtained with the 
two-loop formula with the same A^s value and with the same NLO photon structure functions 
from GRV. In a genuine LO calculation with one-loop a s and LO structure functions the cross 
section is somewhat larger. That the NLO corrections near R = 1 are not very large is due 
to the dominance of the direct contribution which is not changed by the NLO corrections very 
much H || in contrast to the NLO corrections for the double-resolved cross section 

In Fig. 2, 3, and 4 the inclusive two-jet cross sections — plotted, cis ci functions 

of pt, integrated over the rapidities rji and rj2 as in the experimental analysis of the TOPAZ 
H], AMY m and OPAL || collaborations. The conditions for the TOPAZ experiment are 
#max = 3.2°, 21,22 < 0.75 and |?7i|, |t7 2 | < 0.7. The cone radius R — 1 is used in all three 
figures. The comparison with the TOPAZ data |3J] is seen in Fig. 2. The full curve is the sum 
of the direct, single-resolved and double-resolved contributions. The separate cross sections 
are also plotted. From these curves one can see that at the TRISTAN energy the resolved part 
is significant only at very small px- At large px — 8 GeV the direct contribution is dominant. 
But both, the double-resolved and even more the single-resolved cross section make a non- 
negligible contribution. The agreement of the theoretical prediction with the TOPAZ data is 
very good. It is clear from this comparison that the data cannot be explained by the direct 
contribution alone even at the larger px- In Fig. 3 the same curves for the AMY experiment 
are shown. In this experiment the cuts are # max = 13° and 1 77^ | , \ f]2\ < 1. Z\ and z 2 are 
not restricted further. We show again the three contributions, the direct, single-resolved and 
double-resolved cross section. The data points must be compared with the full curve. The 
agreement with the theoretical curve is again quite good. The AMY cross section is larger 
than the TOPAZ cross section, essentially since # max is larger. In the fall of 1995 LEP was 
run at center-of-mass energies of 130 and 136 GeV . During this rather short run the OPAL 
collaboration was able to measure inclusive jet cross sections and analyze them with the cone 
jet finding algorithm ||. The cone radius was also R = 1. In Fig. 4 we show their data 
of the inclusive two-jet cross section together with our prediction (full curve) and the cross 
sections for the direct, single-resolved and double-resolved part. The kinematical constraints 
are 6 max = 1.43°, |r/i|, \t] 2 \ < 1. Because of the larger cm. energy we have taken Nf = 5 and 
the corresponding A^g = 0.130 GeV. At small px the data are very accurate. The combined 
statistical and systematic errors are shown. The agreement of the theory with the data is good 
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over the whole pt range. We see that the double-resolved contribution is much more important 
at this higher energy. It gives the dominant contribution for px < 5 GeV. Even at pt = 16 GeV 
the direct contribution is only 60% of the total. Except for the last high p? point which has a 
larger error the pure direct curve is not consistent with the data. This shows that the other two 
contributions, the single-resolved and double-resolved, are needed to explain the data. The 
OPAL collaboration measured also the inclusive single-jet cross section. The comparison of 
their data with our predictions is shown in M and M. The agreement between theory and the 
data is as good as in the two-jet case. The inclusive two-jet cross section for LEP1.5 contains an 
extra factor 2 as compared to the definition used for the TRISTAN energy in order to account 
for the definition as used in ||. 

In all the comparisons between experiment and theory shown above, we employed the GRV 
presentation of the parton distributions of the photon. Other NLO parton distributions are 
those of Gordon and Storrow [11] and Aurenche et al. [|l^]. Before we compare the data with 



predictions based on these other parton distributions we investigate the scale dependence of the 
inclusive two-jet cross sections for the two e + e~ cm. energies corresponding to TRISTAN and 
LEP1.5. The results are shown in Fig. 5a,b where we have plotted the two-jet cross sections 
for three different scales /i = M = S,Pt with £ = 1/2, 1 and 2, and compared them with the 
TOPAZ (Fig. 5a) and OPAL data (Fig. 5b). In both cases the cross section is remarkably 
independent of the scale change. Only at rather small pt < 4 GeV the cross section in Fig. 5a 
exhibits some scale change. For px > 4 GeV the cross section at = 58 GeV deviates less 
than ±4% from the standard with £ = 1. Similarly the cross section at = 133 GeV shown 
in Fig. 5b varies by less than ±2% above pt = 8 GeV . This small scale change of the 77 jet 
cross sections (the inclusive one-jet cross section shows a similar behaviour) results from the 
fact that the scale change originating from the change of the renormalization scale fi leading to 
a decrease of the cross section with increasing £ is to a large extend cancelled by the dependence 
on the factorization scale M, which is opposite to the change in /1. The dependence on M is 
mostly cancelled between the NLO direct and the LO single-resolved cross section and similarly 
between NLO single-resolved and LO double- resolved cross section J7|. But some dependence 
on M remains when all three cross sections are evaluated in NLO. Thus this reduced scale 
dependence of the 77 jet cross sections would be a good place to determine Ams- For this we 
need photon structure functions with varying A^g which are not available yet and accurate data 
in particular for large px- In order to see the change of the two-jet cross section when we change 
from Ag = 0.130 GeV to A§ = 0.250 GeV" in the formula for a s we have plotted the cross 
section for this large value of A^g- also in Fig. 5b. As to be expected the cross section is larger 
now. At the smaller px < 5 GeV the data points lie below the theoretical curve. Since A^g was 
not changed in the photon structure function consistently we do not draw any conclusions from 
this observation. At large pt the data points are still consistent with the theoretical curve for 
the larger A^-. The equivalent result at \/~S = 58 GeV is shown in Fig. 5a, where we compare 
the results for A® = 0.200 GeV with Aj§ = 0.358 GeV which is equivalent to A§ = 0.250 GeV 
when we choose four active flavours. Due to the larger experimental errors the curve with the 
larger A-^s value is still compatible with the data. At the larger px the agreement is even 
somewhat improved. This shows that with the limited experimental accuracy of the data at 
both energies it seems difficult to draw clear conclusions even if the photon structure functions 
were given for varying A^g. 

In Fig. 6a,b we investigated the influence of the gluon distribution of the photon on the two- 
jet cross section. It is clear that the gluon distribution will have a remarkable influence only at 
small pt since this distribution function is dominant for small x. To see this we have assumed 
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Fgh = with the assumption that the quark distribution functions are already well enough 
constrained by the deep inelastic e— 7 scattering data. As is seen in Fig. 6a,b the prediction 
with Fg/^ = gives very bad agreement with the data at small p T (p T < 5 GeV in Fig. 6a and 
Pt < 8 GeV in Fig. 6b). From this we conclude that the two-jet data require a non- vanishing 
Fg/ry in the photon structure function. This agrees with earlier findings in j|, [5] on the basis of 
more model dependent calculations and with LO calculations in [14] based on new GS structure 
functions. 

The last point is the prediction for other parton distributions of the photon than the GRV 
distribution in NLO. The comparison is shown in Fig. 7a,b together with the TOPAZ and OPAL 
experimental data. As we can see there is little difference between the results for the GRV 



10|| and the ACFGP [JT2J structure functions. The cross sections for the GS |TT[ distribution 
functions are somewhat smaller than for GRV. In average, for both cm. energies, the data 
agree better with the GRV and ACFGP curve than with the curve based on the GS structure 
function, in particular in the comparison with the TOPAZ data. In order to exclude one of 
these structure functions much more accurate data are needed. In general, we conclude, that 
all three NLO parton distributions of the photon account remarkably well for the existing data 
of the inclusive two-jet production. It is clear, that additional information, as for example, the 
rapidity distributions for the two jets could constrain these structure functions further. 

The same comparison has been done with the AMY two-jet data. Since in this case, the 
GRV curve lies at small px above the data points the agreement with the prediction with the 
GS structure function is better than for the TOPAZ data. 

Overall we can state that the agreement between the theoretical predictions and the mea- 
surements of the two-jet cross sections is remarkably good. This is true for the magnitude 
of the cross section and its shape as a function of px- We should not conceal, however, that 
the NLO calculation gives the jet cross sections for massless partons, whereas the experimental 
jet cross sections are measured from hadrons. Assuming the usual parton-hadron duality, no 
difference between parton and hadron level would be observed. This parton-hadron duality 
is very well obeyed in jet production for e + e~ annihilation at the Z resonance. In this case 
the hard scale is larger than in the 77 processes, so that it is problematic to draw any firm 
conclusions from it. On the other hand it is not easy to estimate the hadronization effects for 
jet production in 77 processes. For this one needs a Monte Carlo routine which generates the 
final state partons including the NLO correction with subsequent fragmentation into hadron 
jets. Such a routine is not available yet. Estimates on the basis of LO Monte Carlos may not 
be too conclusive. Further details for this case are reported in M. 



4 Summary 

Differential inclusive two-jet cross sections da/dpT have been calculated as a function of px in 
NLO for the direct, single-resolved and double-resolved contributions. The superposition of 
these cross sections are compared to the TOPAZ, AMY and OPAL data. We obtained good 
agreement between measured data and the theoretical predictions. The same is true for the 
single inclusive cross sections [EL [7| . The single- and double-resolved cross sections are obtained 
with the GRV [|K| photon structure function. Two other sets GS (Tl|] and ACFGP ||T2 | 



were 



also investigated. With them the data are described almost equally well. The 77 jet cross 
sections are remarkably stable against scale changes. With increasing center-of-mass energy 
the double-resolved contribution becomes more and more important. This is clearly seen by 
comparing the prediction for the TRISTAN and LEP1.5 energies. Similar to earlier findings 
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we exclude a description of the data without a gluon distribution in the photon. It is hoped 
that more accurate data from LEP2 may be able to constrain the parton distributions of the 
photon further. 
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6 Figure Captions 



Fig. 1: The inclusive two-jet cross section as a function of the cone radius R for LEP2 kinematics 
and pt = 10 GeV, T]i = r] 2 = 0. 

Fig. 2: Inclusive two-jet cross section da/dpr with R = 1 as a function of pt for direct (dashed), 
single-resolved (dashed-dotted) and double-resolved (dotted) production in NLO. The 
full curve is the sum of all three components compared to TOPAZ data (\rji\, \r] 2 \ < 0.7) 

Fig. 3: Inclusive two-jet cross section da/dpr with R = 1 as a function of pr for direct (dashed), 
single-resolved (dashed-dotted) and double-resolved (dotted) production in NLO. The 
full curve is the sum of all three components compared to AMY data (|^i|, \t]2\ < 1) 

Fig. 4: Inclusive two-jet cross section da/dpr with R = 1 as a function of pt for direct (dashed), 
single-resolved (dashed-dotted) and double-resolved (dotted) production in NLO. The 
full curve is the sum of all three components compared to OPAL data (|^i|, 1^1 < !)• 

Fig. 5: Inclusive two-jet cross sections da/dpT for three different scales fi = M = t^p? with £ = ~ 
(dashed), £ = 1 (full), £ = 2 (dotted) and different A^g- (dashed-dotted) compared to (a) 
TOPAZ data (A^ = 0.358 GeV) and (b) OPAL data (A§ = 0.250 GeV). 

Fig. 6: Inclusive two-jet cross section da /dpx with GRV parton distributions and with F g ^ = 
(dashed) and F 9 / 7 ^ (full) compared to (a) TOPAZ data and (b) OPAL data. 

Fig. 7: Inclusive two-jet cross section da/dpx for different parton distribution sets: GRV [|TIJ 
(full), GS (dashed) and ACFGP (dotted) compared to (a) TOPAZ and (b) 
OPAL two-jet data. 
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